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Abstract 

Thermal leptogenesis in the canonical seesaw model in supersymmetry suffers from the incom- 
patibility of a generic lower bound on the mass scale of the lightest right-handed neutrino and 
the upper bound on the reheating temperature of the Universe after inflation. This is resolved by 
adding an extra singlet superfield, with a discrete Z2 symmetry, to the NMSSM (Next to Minimal 
Supersymmetric Standard Model). This generic mechanism is applicable to any supersymmetric 
model for lowering the scale of leptogenesis. 
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I. INTRODUCTION 



In the minimal standard model (SM), neutrinos are massless. However, small nonzero neu- 
trino masses are required by the atmospheric and solar neutrino experiments. A natural 
explanation for such tiny neutrino masses in the SM comes from an effective dimension-5 

a 



operator 



(1) 



where [va-, la), Oi = e, fi, t are the usual left-handed lepton doublets transforming as (2, —1/2) 
under the standard electroweak SU{2)l x U{1)y gauge group and (0^,0°) ~ (2, 1/2) is the 



usual Higgs doublet of the SM. There are three realizations of this operator 



popular one being the canonical seesaw 



2|, the most 



3| mechanism which adds three singlet heavy neutral 



fermions iVj, i = 1, 2, 3 to the SM Lagrangian. The neutrino mass matrix is then given by 



(2) 



where AAdis the 3x3 Dirac mass matrix linking with Ni through the Yukawa interactions 
The Majorana masses of Ni violate lepton number by two units. Therefore, in the early 



Universe, a net lepton asymmetry may be generated 



J| through the out-of-equihbrium decay 



of the lightest Ni (call it A^^i). The generated lepton asymmetry then gets converted into a 
baryon asymmetry through the interactions of the SM sphalerons 5(1 which conserve B — L, 
but violate B + L, where B and L are baryon and lepton number respectively. The existence 
of Ni explains thus at the same time why both neutrino masses as well as the observed 
baryon asymmetry of the Universe (BAU) are nonzero and small. 

In supersymmetric theories the reheating temperature (T/j) following inflation is likely to 



be rather low 



a, LZt]. Although some models 



81] may allow a higher reheating temperature, in 
the conventional models T/j is bounded strongly from above by the possible overproduction 
of gravitinos. On the other hand, in the simplest version of the seesaw mechanism the 
condition for thermal leptogenesis requires that the mass Mi of the lightest right-handed 



2 



neutrino A^i should be much higher than Th, assuming that A^i contributes to the left-handed 
neutrino masses dominantly. Since inflation would erase any pre-existing lepton asymmetry, 
the asymmetry generated by A^^i after inflation would be highly suppressed by its small 
number density and hence this mechanism will fail to explain the BAU. 

To avoid this problem, several ideas have been discussed in the literature 



10 



11 



13 



121. 



14l |. In this paper we follow 



An attractive scenario is the extended seesaw mechanism 
the same scheme and work with the canonical seesaw mechanism (SM plus three Ni) in the 
Next to Minimal Supersymmetric Standard Model (NMSSM), which has an extra singlet 
superfleld x- To distinguish Ni from x, an exactly conserved Z2 discrete symmetry is im- 
posed, corresponding to (—1)^. We then add an extra heavy singlet superfleld S, together 
with a softly broken discrete symmetry Z2, under which S is odd and all others are even. 
As a result, the production and decay channels of S are different, and the out-of-equilibrium 
decay of S can take place much below the mass scale of the lightest right-handed neutrino 
A^^i. The lower bound on the mass scale of 5* can then be compatible with the upper bound 
on the reheating temperature after inflation. 

The rest of this paper is arranged as follows. In section II we review the canonical 
leptogenesis and briefly recall the Davidson- Ibarra (DI) bound on the mass scale of A^^i. In 
section III we discuss an extended seesaw model by introducing an additional heavy singlet 
fermion S of mass less than that of A^^i. In section IV we discuss how the thermal-leptogenesis 
constraint on the mass scale of 5* can be lowered in comparison with the mass scale of Ni. 
In section V we solve the required Boltzmann equations numerically and show how the 
low mass scale of S is compatible with thermal leptogenesis. In section VI we state our 
conclusions. 



II. CANONICAL LEPTOGENESIS AND DI BOUND 

In canonical leptogenesis the lightest right-handed neutrino Ni decays into either i~(f)'^ and 
z/0°, or i~^(f)~ and Thus a CP asymmetry can be established from the interference of 



3 



the tree-level amplitudes with the one-loop vertex 
lepton asymmetry arises when the decay rate 



~\ and self-energy corrections 15^ 



15 . A net 



= ^^1^1 (3) 
fails to compete with the expansion rate of the Universe 

H{T) = IMgl/'^ (4) 

at T ~ Ml, where rhi = {171^^1710)11/ Mi is the effective neutrino mass parameter, (7* ~ 228 
is the effective number of relativistic degrees of freedom in the MSSM and Mpi = 1.2 x 10^^ 
GeV. This means that a upper bound on rhi may be established by first considering the 
out-of-equilibrium condition H{T = Mi) > Ti which gives 

mi < 1.6 X 10"^e1/. (5) 

However, for rhi > lO^^eV a reduced lepton asymmetry may still be generated, depending 
on the details of the Boltzmann equations which quantify the deviation from equilibrium of 
the process in question. 

Assuming a normal mass hierarchy in the right handed neutrino mass spectrum the CP 
asymmetry ei is given by 

^ ^ 3_ /MA Jm[(ml,m^)i2]^ 

' Sttv^ \M2j (mj^mz?) 11 ' 



16l ] with precision. 



The baryon-to-photon ratio of number densities has been measured 

i.e. 

r]B-^ = 6.ll°ixlO-^°. (7) 



To get the correct value of rjB, one needs ei to be o 
DI upper bound on the CP asymmetry parameter 



' order 10 ^ to 10 ^. However, using the 
9| 



, , 3Mi r- ^ 

|ei| < ^^V^^^atm (8) 

one can get a lower bound on the mass scale of the lightest right-handed neutrino as 

q^../ Vb \/4x10-3\/ V \2/ 0.05eV \ 



where v is the vacuum expectation value (vev) of the SM Higgs and it is also assumed that 
the physical left handed neutrinos follow the normal mass hierarchy. Since is not likely 
to exceed 10^ GeV, this poses a problem for canonical leptogenesis. In order to overcome 
this problem we consider an extended seesaw model as follows. 

III. THE MODEL FOR THERMAL LEPTOGENESIS BELOW THE DI BOUND 



In a recent paper 13| we proposed a singlet mechanism to overcome the DI bound shown 
in Eq. ([9]). We now present a realistic model where this mechanism can be implemented. 
One more ingredient has been added to produce these singlet fields abundantly in a thermal 
bath. 

We start with the NMSSM model, which includes a singlet superfield x addition to 
the usual particles of the Minimal Supersymmetric Standard Model (MSSM). To implement 
the seesaw mechanism, we also include three right-handed neutrinos A^j, i = 1,2,3. We 
then demonstrate how a minimal extension of this model may admit a very low scale of 
leptogenesis, thus overcoming the gravitino problem. We include another singlet superfield 
S and impose a x discrete symmetry. Under Z2, the lepton superfields Lj, l^, Ni, S are 
odd, whereas the Higgs superfields 0i^2, X ^-^'g even. This corresponds to having an exactly 
conserved lepton number (—1)^, or the usual i?— parity of the MSSM. Under Z'2, S is odd and 
all others are even, but Z2 is allowed to be broken softly. The most general superpotential 
invariant under Z2 x Z2 is then given by 

W = hljLirj(t>i + hi,UNj(t>2 + ii<pi(t>2 + Mi.NiN, + M^XX 

+OiXXX + + fNXN^N^ + MsSS + fsxSS. (10) 

We do not discuss quarks nor other interactions of the NMSSM, which have been studied 
elsewhere. We deal only with neutrino masses and leptogenesis. 
We now break Z2 softly and the only possible such term is 



5 



i.e. exactly as required to implement the singlet mechanism of ref. 13]. This allows 5* to 
mix with iVj to form a 7 x 7 mass matrix in the basis [Li Ni S], i.e. 



mo ^ 



M 



mo Mm d 
d Ms 



(11) 



where Ms = fsix)^ without loss of generality we choose M^ to be diagonal with masses 
^1,2,3- For small di/{Mi — Ms) as well as the usual seesaw assumption that the entries of 
m£, are very small relative to Mj, the heavy masses are roughly given by 



Ms' 



M, 



Mi + 



M,-Ms'' 
d^ 



Mi -Ms' 

corresponding to the mass eigenstates 5" and 

di 



(12) 



5' ~ 5 - 5^ 



iv; ~ AT. + 



Mi -Ms 



-Ni 



di 



M,-Ms 



-S. 



(13) 



The light neutrino mass matrix is then 



[m„)ij ^ 



dl 



M. - M, 



{mo) 



kj- 



(14) 



In the limit di ^ we recover the neutrino masses as in the canonical seesaw mechanism. 
We assume thus Ms' — Ms and M^' ~ Mn in the following. 



IV. LEPTON ASYMMETRY AND LOWER BOUND ON Ms 

In this model, the addition of S allows the choice Ms < Mi. The induced couplings of S 
to leptons are suppressed by factors of di/Mi compared to those of A^,. The decay rate of 5* 
is thus given by 

rf, = ^Ms ^ [mM,{d,/Mif] , (15) 



where we have assumed Ms -C Mi and the effective neutrino mass parameter is defined as 

(m|)mD)n 



Mi 



(16) 



Assuming that -^l ~ ~ ^"^^ ^3 ^ ""^2 > the above Eq. f|T5l) can be rewritten as 



D — 



-Msm^M^{d^/M^f . 



(17) 



The out-of-equihbrium condition Ff, < if(T ~ Ms) is thus suppressed by a factor 



(i| \ / _ / m3 



(18) 



in comparison to Ti/H{T ~ Mi) and hence can be satisfied at a lower mass depending on 
the value of d^. The value of ms, fhi and M3 can be approximately fixed from the low energy 
neutrino oscillation data. So, the remaining free parameters are cis and Ms on which the 
suppression factor rj depends. 



S N 1 N 



N 



S N 1 




FIG. 1: Tree-level and one- loop (self-energy and vertex) diagrams for 5 decay, which interfere to 
generate a lepton asymmetry. 



The CP asymmetry generated by the decays of S comes from the interference of the 
tree-level and one- loop diagrams of fig. ([1]). Both the numerator and denominator of Eq. 
([6]) are then suppressed by the same {di/MiY factor, and we obtain 

3 (Ms\ Im[{m^j,mD)i2f 



(19) 



11 



where we have assumed that Mi <^ M2 <^ M3. Thus we see that the CP asymmetry 
parameter is independent of the suppression parameter ^3. Therefore, depending on the 
value of K the L-asymmetry will saturate at different temperatures as implied by Eq. f|T8l) . 
This is shown in Section V by numerically solving the required Boltzmann equations. As 
demonstrated in fig. ([H]) the value of Ms can be lowered much less than the DI bound on 
Ml by an appropriate choice of d^. This is because the low values of Ms are not restricted 
by the low energy neutrino oscillation data for di -C Mj as we have seen from Eq. (12). 
Moreover the washout effects are suppressed for low values of {d^/M^). So, a successful 
lepton asymmetric universe before the electroweak phase transition can be created even for 
a TeV scale of S. 



V. NUMERICAL ESTIMATION OF LEPTON ASYMMETRY 
A. Production and decay of S 

In this model S is produced through the decay of x- The corresponding Yukawa coupling 
fs can be as large as of order unity. Therefore, 5* can be brought to thermal equilibrium 
through the scattering processes: SS SS, SS XX^ ^'^'^ xS xS- Note that these 
processes never change the number density of S, but they keep S in kinetic equilibrium. 
The decay rate of x can be given by 

where {K1/K2) is the boost factor and z = Ms/T. Thus the inverse decay of x is given by 

Once the S particles are produced, they decay through the channel: S i(f)^ , as shown 
in fig. (JH) which violates lepton number by two units. Apart from that, the other process 
which depletes the number density of S is Si (f> Qt. This is shown in fig. ([2]). 

The subsequent decay of 5*, below its mass scale, then produces the required baryon 
asymmetry through the leptogenesis route. However, an exact lepton asymmetry can be 



8 



FIG. 2: AL = ±1 processes which deplete the number density of S. These processes also deplete 
the net lepton number density produced through the decay channel. 



17l |. It is useful to define the Boltz- 



estimated by solving the required Boltzmann equations 
mann equations in terms of the dimensionless variables Ys = ns/s and Yl = til/s, where 
Ys is called the co moving density of S while Yl is the density of net lepton in a comoving 
volume and 



45 



(21) 



is the entropy density. The required Boltzmann equations are given as 



and 



dYs 
dz 



dYr 



^ D 



+ 



_zH{z) zH{z)_ 



(-S- 



(Ys - Yl'') - 



(22) 



(23) 



dz " zH{z) ' ' zH{z) 

where r|,, Ff and F^ simultaneously represent the decay, scattering and wash out rates 
that takepart in establishing a net lepton asymmetry in a thermal plasma. The Hubble 
parameter H{z) is given by 



Hiz) 



with H{Ms) = l-Q7gl/^j/- 



(24) 



In a relativistic frame the decay rate f[T7l) can be rewritten as 



Ko(z) 



(25) 



The Ff in Eq. (!22|) represents the processes which violate lepton number by one unit and 
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FIG. 3: AL = ±2 processes which deplete the number density of net leptons. 



is given by ^ 

Ts = 4rJ, + 2T' . (26) 



The T]v in Eq. represents the lepton number violating processes by two units and 
is given by 

= 2^/1) + 2r|^,t + 2r^ (^y^^ + 2r^i + 2r^i t > (27) 

where Ff^ = (n^^/n^^)F£,. In Eqs. ( I26l) and ( 127|) the F's are defined as F^ = ipfi/rfx) where 
7j is the scattering density. Note that the other AL = ±2 processes: II — > NiSNi —>■ (jxf), 
and of course higher order processes, which contribute to F^ are suppressed in comparison 
to the processes shown in fig. ([3]). 

B. Solution of Boltzmann equations 

In fig. (jl]) we have plotted the decay and inverse decay of x and S against z. It is shown that 
the inverse decay of £ + 0^ ^ is not sufficient to bring S into thermal equilibrium even 
if the suppression factor ^3 is as large as 10® GeV. On the contrary, the decay rate of x is 
sufficiently larger than the Hubble expansion parameter. Hence S can be brought to thermal 
equilibrium through the scattering process involving x as long as (M^/Ms) — 0(10^^^). 
Therefore, at a temperature far above the mass scale of S it is in thermal equilibrium and 
hence a net lepton asymmetry in the thermal plasma is zero. Below the mass scale of 5* 

^ We have not included the SUSY processes. It is shown that upon inclusion of those processes the result 
doesn't change significantly [3]. 
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FIG. 4: The production and decay, and their inverse processes, are compared with the Hubble 
expansion parameter H for a typical set of parameters. We have used Ms = 10* GeV, = 10^*^ 
GeV, M3 = m^^GeV, = lO^^ eV, = 10* GeV and fs = 0.5. 

the lepton number violating processes go out of thermal equilibrium and thus produce a net 
lepton asymmetry dynamically. This is obtained by solving the Boltzmann equations (l22l) 
and (123|) . We take the following initial conditions: 



The evolution of the number density of 5* and the corresponding asymmetry with respect 
to z are shown in figs. ([5]), ([6]) and ([7]). At any epoch z the value of Ys and the corresponding 
asymmetry Yl can be inferred from 



Since the decay rate of S is suppressed by a factor of the asymmetry is produced at late 
times depending how small it is. However, the value of k cannot be made indefinitely small. 
Because a net L-asymmetry has to build up before the electroweak phase transition which 
is required to be converted to the 5-asymmetry through the sphaleron transitions. The 



Ys = Yl'' and 1^ = at z^O. 



(28) 




z 



(29) 
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FIG. 5: The evolution of S is shown against z with Ms = 10^ GeV, Mi = 10^ GeV, M3 = 10^° 
GeV and ^3 = 10^ GeV and the CP asymmetry parameter is = 10~^. 

final L-asymmetry is numerically obtained for three values of ^3 in figs. ([5]), ([6]) and ([7j). 
It is found that the final L-asymmetry is almost same apart from a numerical factor. This 
is because for the delayed decay of 5* the wash out effects are comparatively small. While 
K = 10~^ and 10""^ are used in figs. ([5]) and ([6]), it is of 10"^ in fig. ([7]). As seen from figs. (Q, 
dED and (C]), for k = 10"^ 10-^ lO'^ the value of Yl is saturated at around Zs ~ 10, 10^ 10^ 
respectively. Assuming that a final L-asymmetry has to be produced before T^w ~ 100 GeV 
the minimum tolerable value of k = 10^^^ is obtained. This indicates that for M3 = 10^° 
GeV, (ig and Ms can be readjusted among themselves so as to get the suppression factor k, 
ranging from 10~^ to 10^^^. This is shown in fig.®. The solid line in fig. ([8]) is obtained 
for Ms = ds- The region above to that are defined by ^3 > Ms- So these values of ^3 
are unnatural and are not allowed. While the region below to the solid line are defined by 
< Ms and hence is allowed for naturalness. Thus we see that a wide range of Ms values 
from 10'^ GeV to 10*^ GeV are allowed that can produce the required lepton asymmetry 
before the electroweak phase transition. 
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FIG. 6: The evolution of S is shown against z with Ms = 10^ GeV, Mi = 10^ GeV, M3 = 10^° 
GeV and ^3 = 10'' GeV and the CP asymmetry parameter is = 10~^. 

VI. CONCLUSIONS 

We accomplished the baryogenesis via leptogenesis from the decay of an additional singlet 
5* in a supersymmetric extended NMSSM. The bound coming from the out-of-equilibrium 
condition could be evaded because the couplings of the singlets cancel out from the asym- 
metry, so the couplings could be small and can satisfy the out-of-equilibrium condition even 
at low scales. In the simplest seesaw models the couplings of the lightest right-handed 
neutrino could not be lowered much because that will not enable the thermal production 
of these fields. However, in the present case there is one additional singlet field (x) which 
can produce these 5* fields, having large couplings to them, but itself not taking part in 
leptogenesis. 
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FIG. 7: The evolution of S is shown against z with Ms = 10^ GeV, Mi = 10^ GeV, M3 = 10^° 
GeV and ^3 = 10^ GeV and the CP asymmetry parameter is = 10~^. 

APPENDIX A: SCATTERING DENSITIES 



In this appendix we give the various scattering densities that have been used in the 
Boltzmann Eqs. ([22]) and ([23]). 



70,. 



70,i 



1 

xi xi \ y 



(Al) 
(A2) 



where v is the vev of SM Higgs and Xi = s being the Mandelstam variable, and y 



INI = / dX2y/X2Ki Zy/X2- 



Mr 



'Ms 



1 



X2 



lNl,t 



D,{x2) ' 2Dl{x2) 
MfMsmj 



dX2\fx>2Kx z^fx>2 



Ml 



^i^ + lln(x2 + l) 

X2 + I X2 



where X2 = ^ and 



X2-I 



(A3) 
(A4) 

(A5) 
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FIG. 8: The allowed values of Mg are shown in the plane of (i| versus M5. 
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